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Abstract: For SU(2) gauge theory on the three-sphere we implement the 
influence of the boundary of the fundamental domain, and in particular the 
^-dependence, on a subspace of low-energy modes of the gauge field. We 
construct a basis of functions that respect these boundary conditions and 
use these in a variational approximation of the spectrum of the lowest order 
effective hamiltonian. 
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in : 1 Introduction 



It is our aim to study the dynamics of the low-energy modes of pure SU(2) gauge theory 
defined in a finite volume []T|, ^ ^ . If we take the volume small, the mechanism of asymptotic 
freedom results in a small coupling constant and we can use standard perturbation theory. 
Gradually increasing the volume then allows us to monitor the onset of non-perturbative 
phenomena. Especially when the non-perturbative effects manifest themselves appreciably 
only in a small number of low-lying energy modes, this can be described adequately using 
a hamiltonian formulation. We are interested in the influence of the multiple vacuum 
structure of the theory on the glueball spectrum: in particular we would like to see the 
dependence of the energies on the #-angle which corresponds to the transition over instanton 
barriers. 

In our approach, we impose the Coulomb gauge by restricting the gauge fields to a 
so-called fundamental domain [|], [5| . This is a convex subset of the space of all gauge field 
configurations that is in one-to-one correspondence (modulo constant gauge transforma- 
tions) with the space of gauge orbits. The latter is precisely the physical configuration 
space on which we want to study the dynamics. 

Our strategy will be to split up the gauge field in orthogonal modes and to reduce the 
dynamics of this infinite number of degrees of freedom to a quantum mechanical problem 
with a finite number of modes. We do this by taking standard harmonic oscillator wave 
functions for the high-energy modes of the field. For the low-energy modes, the onset of non- 
perturbative effects means that the wave functional starts to spread out over configuration 
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space and will, in particular, become sensitive to conditions imposed on the boundary of 
the fundamental domain. We therefore replace the wave functionals for these modes by 
functions that satisfy the (^-dependent) boundary conditions. 

The dynamics of these modes can be regarded as an effective low-energy theory. Similar 
to a Born-Oppenheimer approximation we integrate out the fast (high-energy) modes and 
are left with an effective theory of the slow modes. If we are at energies at which the higher 
modes cannot be excited, these modes will, through virtual processes, in first order only 
result in a renormalisation of the coupling constant. 

To make contact with lattice calculations it would be most natural to take the finite 
(spatial) volume to be a 3-torus T 3 . Detailed knowledge of the vacuum structure in this 
case is however limited. In particular the instantons, which are the gauge field configu- 
rations that describe tunnelling between different vacua, are only known numerically ||. 
To circumvent this problem we take our space to be the three-sphere S 3 on which the 
instantons are known analytically. 

Previous studies of gauge theory on S 3 M did not use conditions at the boundary of 
the fundamental domain, but focussed on the influence of the Gribov horizon: rescaling 
the wave function with the square root of the Faddeev-Popov determinant results in a 
potential with an infinite barrier at the Gribov horizon. 

In this paper, we first construct suitable boundary conditions on the space of low-energy 
modes and then give a basis of functions that respect these boundary conditions. Using 
this basis we perform a Rayleigh-Ritz analysis to approximate the spectrum of the lowest 
order effective hamiltonian. The obtained energy differences give a first hint at the values 
of the glueball masses in the different (scalar, tensor) sectors. 

2 The Effective Theory 

We will briefly review the formalism developed in H, ||| for gauge fields on S 3 and construct 
the effective hamiltonian. We embed S 3 in IR 4 by considering the unit sphere parametrized 
by a unit vector n M . Using the scale-invariance of the classical hamiltonian, we can make 
the restriction to a sphere of radius R = 1. -R-dependence can be reinstated on dimensional 
grounds. We introduce the unit quaternions <r M and their conjugates = cr* by 

cr M = {id, it), = (id, -if). (1) 

They satisfy the multiplication rules 

°\&v = V^a, = V^a, (2) 

where we used the 't Hooft r\ symbols ||, generalised slightly to include a component 
symmetric in /i and v for a = 0. We can use r] and fj to define orthonormal framings of S 3 , 
which were motivated by the particularly simple form of the instanton vector potentials 
in these framings. The framing for S 3 is obtained from the framing of 1R 4 by restricting 
in the following equation the four-index a to a three-index a (for a = one obtains the 
normal on S 3 ): 

^ = rj° v n u , e° = f]°„n v . (3) 
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A gauge field on S (n^A^ = 0) can be written with respect to either framing, e.g. 

A, = M = Ay^ (4) 

A gauge transformation g acts as follows: 

{'A\ = g-'Ag + g~ l d ig = g- l D\ A) g. (5) 



In order to isolate the lowest energy levels, we define the quadratic fluctuation operator 
M by 

V(A) =~J sZ \ tr (4) = tr (AiMijAj) + O (A 3 ) . (6) 

The 18 dimensional space A(c, d) given by 

Mc,d) = (c^ + d^)^ (7) 

is the eigenspace of Ai corresponding to its lowest eigenvalue 4, whereas the next eigenvalue 
is 9. As was explained in @, [J, this space contains gauge copies of the vacuum A = 
and the tunnelling paths inbetween these vacua. In particular, the vacuum A a = —o a 
(c" = —25", df = 0) is a copy of A = under the gauge transformation g = n ■ a with 
winding number one. The tunnelling path is cf = —uSf, df = with u running from to 2. 
For u — 1 it passes through the sphaleron, which is a saddle point of the energy functional. 
The vacuum c\ = 0, df — —25f is a copy of A = under g = n ■ a. The corresponding 
saddle point is called " anti-sphaleron" and is actually a copy of the sphaleron under this 
gauge transformation. The energy functional for these 18 modes is given by 

V(c, d) = - \ tr (4) = V(c) + V(d) + ^ {(4) 2 (4) 2 - «^i) 2 } > (8) 

V(c) = 2tt 2 {2(c«) 2 + 6detc+ i[(c^) 2 - (c^ S a ) 2 ]} . (9) 
The lowest order hamiltonian for these modes is 

wMfe+^)+r^ (10) 

with / = and g(R) the renormalized coupling constant. 

The boundary of the fundamental domain in the full (c, d) space is not known, but at 
the energies we are interested in, the boundary conditions are relevant only at those points 
where the potential is low, i.e. at the sphalerons. At other boundary points the potential 
will be (much) higher then the energy of the wave functional. This means that the wave 
functional will have decayed exponentially at these points and that the precise boundary 
conditions should not have a big influence on the spectrum. By the same token, the precise 
location of the boundary in these regions is not important either. This gives us the freedom 
to choose tractable boundary conditions. 

At the sphalerons however, the boundary conditions are fixed. Since the gauge trans- 
formation connecting the two sphalerons has winding number one, we have to set 

^(A(Sph, 0)) = e ie V{A(0, Sph)), (11) 

thus introducing the #-angle. 
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3 The Rayleigh-Ritz basis 



We now want to determine the spectrum of the hamiltonian (|10D, properly restricted to 
the fundamental domain. The variational basis we need must incorporate the boundary 
conditions but also respect the symmetries of the hamiltonian to obtain an optimal block 
diagonalisation. 

We first turn to the boundary conditions. We define radial coordinates r c and r<i by 



i i 



r c =[cX?, r d = [d°d°\*. (12) 

The sphaleron has radial coordinates (\/3, 0) and angular coordinates c? = —5" 1 (with 
cf = cf/r c ). It will be connected with the anti-sphaleron at (0, a/3). We will restrict the 
( r c r d) plane by r c < a/3, < a/3 and impose boundary conditions at the edges. This 
means that we will be working towards basis functions of the form 4>(r c , rd)Y(c, d). 

As was explained earlier, for values of the coupling constant at which our approximation 
will be valid, only the effect of the boundary conditions at the sphaleron will be felt. By 
imposing boundary conditions in the (r c , r^) plane we pair up two submanifolds, of which 
only the sphaleron/anti-sphaleron need belong to the boundary of the fundamental domain. 

Consider the following decomposition of the full wave function 

^iz^fe^aw. (is) 

c d n 

where q denotes all the modes orthogonal to the c and d modes. The functions xCai^) 
are chosen to be (perturbative) eigenstates of the transverse hamiltonian. In the adiabatic 
approximation we assume the transverse wave function to be in its ground state and we 
assume this ground state to decouple dynamically from the excited states. This results in 
a hamiltonian for ip = ip^ given by 




Q r 2 Q r 




2 + 7^2 ~ 12 I - + - I + + ~ A * I + ~? (14) 



with the laplacian in the angular coordinates. The boundary condition on ip follows 
directly from (|TTJ) : 

^(Sph,0) =e w V(0,Sph). (15) 

To obtain the condition on the derivative of ip we must take the residual constant gauge 
symmetry into account. Focussing on the sphaleron path c" = —uSf, we take coordinates 
c = S(a)H(u,hi), with S G SO(3) and H symmetric. Under a constant gauge transfor- 
mation only S changes: H(u, hi) forms a set of gauge invariant curvilinear coordinates. 
Matching along the sphaleron path across the boundary of the fundamental domain, we 
need to compensate for the curvature with the appropriate jacobian factor u 3 ^ 2 (details 
will appear elsewhere), which results in the following boundary condition on ip: 

— (Sph,0) = -e — (0,Sph). (16) 

We now turn to the symmetries of the hamiltonian. 7i(c, d) is invariant under the trans- 
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(c\L; l s , Z r , T] m s , m r ) 


(cd\ j, m, l s 


; Li, li, T\\ L 2 , l 2 , T 2) 
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[L S c f 


Uh + 1) 


(L S C ) 2 


Uh + i) 






{L s d ? 


l s {l s + l) 






(J s r 









Ji 







l r {lr + 1) 




Zl(zl + 1) 




m r 


&? 


h{h + i) 






j7 


jU + 1) 

m 




—L(L + 7) 




-ii(ii + 7) 






A," 


-L 2 (L 2 + 7) 



Table 1: Behaviour of the functions of the angular variables 

formation c — > ScRi, d — > SdR 2 with S, R±,R2 G SO (3) and under the interchange c <->• d. 
The generators of left- and right multiplication are Lf, L s c , L d and Lf . The following set 
of operators commutes: 

{H, I s , J R , (L?f + (Lf) 2 , (Lf) 2 + (Lf) 2 ,p} . (17) 

V is defined by Vf(c, d) = f(d, c). On S 3 it corresponds to the parity (n , fi) <->• (no, —ft). 
The operator J s = Lf + Lf implements constant gauge transformations: we have to 
demand (J 5 ) 2 = for physical wave functions. The operator J R = L R + L d is the rotation 
operator. The different sectors under this symmetry correspond to scalar glueballs, vector 
glueballs, etc. 

An orthonormal basis of functions of c is given by the set {(c\L; l s , l r , r; m s , m r )}. Each 
of these functions is a certain polynomial in c, homogeneous of degree L. Its eigenvalues 
under the various symmetries are collected in table |1|. A possible degeneracy is labeled 
by r. Using algebraic manipulation programs, we were able to construct these functions 
explicitely for L < 10. Using Clebsch-Gordan coefficients, we define a function Y llt2 (c, d) 
which is an eigenfunction of J R and of J s (i denotes a set (L, l s , l r , r)): 

Y hi2 (c,d) = (cd\j,m,ix,i 2 ) 

= E E E {-i) h ' h+m M+i ( h h 3 ) 

^ ^ v ) v J \ mi Too -to / 

m s =— l s mj=- ti rri2=— 12 \ / 

(-l)'*~ m *-T=P===f (c|Zu; Z a , Zi, n; m a , m x ){d\L 2 \ l s , Z 2 , r 2 ; -m s , m 2 ). (18) 
y + 1 

Note that in order to have J s = 0, the functions of c and d need to have the same l s . This 
restricts the possible combinations of i\ and i 2 . 

Let us focus on the case 8 = 0. In this case we can in principle take the basis functions 
to be the eigenf unctions of the kinetic part of the hamiltonian. These functions are given 
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by 

^\r c )^(r d )Y™(c,d), (19) 

where L k is determined by z& and <£>jp(r) = 7^J(3+l)(7^) with j p (z) the spherical Bessel 
function of order p. During the variational stage of the calculation, however, the use of 
Bessel functions of different order will lead to a large number of integrals. Therefore we 
take the radial functions to be independent of L\ and L2 and define 

(c, d) = ^ 71 (r c )^ 72 (r d )y ili2 (c, d), (20) 

with v? 7 (r) = 7^3 ( r yr). The functions ^^(c, d) now have singularities at certain points, 
but these are not felt during the variational calculation. Taking even and odd combinations 
gives 

(c, d) = ^ 2 ^ d )+P ^ c ) ( 21 ) 
We implement the boundary conditions ( |1"5|) and ( |1~6| ) for 9 = by imposing the following 
conditions on 71 and 72: 

p=-l : ^(1/3) = Vv (V3) = 0, (22) 

p=1 : ^VJ (V 3 ) = ^V,) (V 3 )=0 . (23) 
Or or 

These conditions are expected to be accurate as long as the wave function transverse to 
the sphaleron path (near the sphalerons) is predominantly in its ground state. For general 
9 we multiply (c, d) with a phase factor exp(i9a(r c , r d )). The function a is a kind 
of Cherns-Simons functional that gives the right behaviour to the wave function under 
large gauge transformations. The resulting functions no longer have well-defined parity, 
but they do obey the general boundary conditions for suitable a. Also the hermiticity of 
the hamiltonian for these functions can be checked explicitely. Sufficient conditions on a 
are: a(r c ,r^) = — a(rd,r c ) and a(\^3, 0) = ±. We choose 




(24) 



For (3 — > 00 we approach the situation that the phase factor over the entire edge is constant 
and equal to e %e . But already for the choice (3 = 2, the boundary conditions at the 
sphalerons are taken into account properly. 



4 The Rayleigh-Ritz analysis 

To perform the Rayleigh-Ritz analysis, we take a finite subset of the basis constructed 
above. We then calculate the matrix of Ti with respect to this truncated basis and di- 
agonalise. The eigenvalues found are upper bounds for the true eigenvalues. In order to 
get lower bounds, we use Temple's inequality || which as input requires the expectation 
values of Ti 2 for the variational eigenvectors. For the calculation of the matrix elements 
we essentially need the radial integrals 

/ dr r n <y9 7 /(r)y3 7 (r) (25) 
Jo 
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with 7 and 7' given by the boundary conditions (|23|) . The angular sector is more difficult: 
especially V 2 contains some non-trivial angular operators. By virtue of the construc- 
tion flUf) all the angular integrations over c and d can be reduced to integrations over just 
c. Both the radial and the angular integrals were tabulated and stored for later use in a 
Fortran program. 

Due to the symmetries, the function space splits up in different sectors characterized 
by the conserved quantum numbers j, m, + 1) + l 2 (h + 1) and, for 9 = 0, p. It is an 
amusing fact that the only influence of the conserved numbers j and m on the spectrum is 
through the number (— l)- 5 and through the fact that /1, 12 and j have to satisfy a triangle 
inequality. If we denote a sector with a triple (/i,/ 2 ,j), this for instance means that the 
sectors (1, 1, 0) and (1, 1, 2) are degenerate. 

We found that for the window of coupling constants for which we trust our approxima- 
tion, the lowest-lying scalar (j = 0) and tensor (j = 2) levels are found in respectively the 
sectors (0, 0, 0) and (1, 1, 2). See figures [l] and @. For 9 = 0, the vacuum corresponds to the 
ground state of the (0, 0, 0)-even sector. The scalar glueball + can be identified with the 
first excited state in the (0, 0, 0)-even sector. Note that for small / the lowest level in the 
(0, 0, 0)-odd sector 0~ has actually a lower energy than the scalar glueball. For small /, 
there is virtually no dependence of the masses on 9, whereas for larger / the 6*-dependence 
is shown in figure H. 



5 Conclusions 

Using the Rayleigh-Ritz method we can determine the spectrum of the effective hamil- 
tonian. The use of Temple's inequality gives us confidence that our results are accurate, 
especially since experience tells us that the actual error is usually much smaller. The results 
are also consistent with an ordinary perturbative calculation of the levels that was done 
up to O (f 3 )- When we make a preliminary comparison with results on the torus, we find 
that the mass ratio m2+/mo+ comes out rather nice: for a volume of 1 fm 3 , corresponding 
roughly to / = 0.4, we obtain the value 1.5. 

Our model is valid for values of the coupling constant at which only the boundary 
conditions at and near the sphalerons are felt. Checking this a posteriori with the help of 
plots of the wave function in the (r c , r^) plane indicates that / should not be larger than 
roughly 0.4. Consider the function 

mr c ,r d )\ 2 = J dcdd\^{c,d)\ 2 (26) 

which is a measure of the probability distribution in the (r c , r<j) plane. Dividing this 
iV'tjo r d) \ by r\r\, we obtain a function with the expected behaviour of the true wave func- 
tion: it is localised at the sphalerons and decays exponentially in the transverse directions. 
Although the lowest barrier is the sphaleron, the volume factor r^r\ causes the configu- 
rations that are close to the sphaleron but have a somewhat higher energy to make the 
dominant contribution to the tunnelling. The relevant parameter here is the characteristic 
decay length of the wave function, which in turn is determined by the rise of the potential 
in the transverse directions. 
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Although we are primarily interested in the case 9 = 0, an appreciable dependence 
on 9 for a certain value of / shows that the non-perturbative influence of the boundary 
has become important. To explain the fact that the spectrum is not exactly periodic in 
9, note that our implementation of the ^-dependence (|24]) only has this periodicity in the 
limit (3 — > oo. The volume effect described above implies that the relevant distribution 
|V , ( r c) r d)| samples a piece of the boundary over which the phase difference already starts 
to depart from e ld . 

The first improvement that has to be made is to take the higher order corrections of 
the other modes into account. The leading order taken here gives rise to a renormalisation 
of the coupling constant. Higher order corrections may alter the potential in our effective 
model and can for instance result in a different degree of localisation around the sphalerons. 
We hope to report on these matters in the near future. 
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Figure 1: Lowest energy levels for 6 = 0. Drawn curves correspond to levels in the (0, 0, 0) 
sector. The dashed curve denotes the ground level in the (1, l,2)-even sector. The short- 
dashed curves are the perturbative expansions, and the individual dots are lower bounds 
on the levels as obtained by Temple's inequality. 
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Figure 2: Glueball masses for 6 = as a function of the coupling constant. The drawn 
curves are the masses of resp. the first scalar (0 + ) and tensor (2 + ) glueball. The dashed 
lines denote the lower bounds, the dotted lines the perturbative result. 
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Figure 3: Scalar glueball mass at / = 0.4 as a function of 9. 
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